Three-quark clusters at finite temperatures and densities 
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We present a relativistic three-body equation to study correlations in a medium of finite temper- 
atures and densities. This equation is derived within a systematic Dyson equation approach and 
includes the dominant medium effects due to Pauli blocking and self energy corrections. Relativity is 
implemented utilizing the light front form. The equation is solved for a zero-range force for param- 
eters close to the confinement-deconfinement transition of QCD. We present correlations between 
two- and three-particle binding energies and calculate the three-body Mott transition. 
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Exploring the phase structure of quantum chromody- 
namics (QCD) for the whole density-temperature plane 
is a challenging task for the standard model. Lattice 
Monte Carlo simulations have revealed several exciting 
results over the past decade In addition, at finite 

densities, modeling of QCD has added substantially to 
our understanding of its rich phase structure. 

One particularly interesting result is the possibil- 
ity of "color-superconductivity" , extensively discussed 
recently The possible appearance of color- 

superconductivity is related to similar ideas that led to 
nucleon or Cooper pairing. Another important aspect 
is the transition from nucleons to quarks as relevant de- 
grees of freedom (Mott transition) that is likely to be 
related to the "confinement-deconfinement" phase tran- 
sition. In these cases the influence of three-quark corre- 
lations has hardly been investigated. Recently, Pepin et 
al. addressed the question of a possible competition of 
three-quark clusters and two-quark condensation at 
finite density but zero temperature. 

The problem of correlations in medium is tackled 
within a Green functions formalism jsj using the Dyson 
equation approach j^] . Three-particle correlations on the 
basis of a new in- medium Alt-Grassberger-Sandhas |ld] 



(AGS) type equation have been investigated in [11 



and implemented into a Boltzmann-Uehling-Uhlenbeck 
simulation of heavy ion collisions and compared to ex- 
periments in [p^l^p^. Recently, the formalism has been 
extended to four-particle correlations and solved to study 
in- medium properties of the a-particle p9| . In contrast 
to the nucleonic phase, however, investigations of the 
quark phase of QCD entail several new issues such as 
relativity, color degrees of freedom, confinement, effective 



residual interactions, number of flavors, among others. 

In order to derive effective few-body equations for 
quarks at finite temperatures and densities we have 
to consider in-medium effects as well as relativity. In 
medium the mass of the constituents may vary and in 
fact becomes close to zero for higher temperatures and 
densities, see e.g. . To implement relativity we use the 
light front approach ||2l| . In the context of the quantum 
statistical framework it has some compelling advantages, 
although the concept of quasi particles in a background 
field that will be utilized further down and extended to 
include correlations introduces a special frame of refer- 
ence. These are: i) Several boosts are kinematical (in- 
teraction free). As a consequence of the transitivity of 
kinematical boosts the Fock state decomposition is sta- 
ble . Since we decouple the hierarchy of Green func- 
tion equations using a cluster mean field expansion 
this decoupling of Fock spaces is therefore retained, u) 
The dependence of the equations for the isolated system 
on the cm. momentum of the cluster can easily be sep- 
arated, see e.g. p3|-p5|. In a homogeneous medium this 
is also an important feature of the light front form, be- 
cause inclusion of Pauli blocking factors that depend on 
the cm. momentum of the cluster leads to a parametric 
dependence on the cm. momentum only, iii) Pair cre- 
ation processes are likely to be suppressed on the light 
front p^,^ in some frames of reference. Therefore we 
may presently consider particles only and no antiparticles 
which leads to technical simplifications that, however, 
can be relaxed as the investigation goes on. iv) Another 
advantage is that the light front form can be formulated 
in a Hamiltonian language and therefore all the results 
of the previously tested quantum statistical framework 
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used here can be recovered using a proper interpretation. 

On the other hand the vacuum structure in the hght 
front dynamics involves technical difficulties in the pres- 
ence of condensates, viz. zero modes in general, see 
e.g. These difficulties include a proper description 

of zero modes and Goldstone modes which are being in- 
vestigated by several groups, see e.g. Further, ro- 
tational invariance on the light front implies interaction 
dependent generators. As a consequence angular decom- 
position is technically more involved. Presently, we make 
use of angular averaging. This introduces uncertainties, 
but for homogeneous infinite matter angular averaging 
has proven quite useful in the past. Finally, the descrip- 
tion of spin degrees of freedom on the light-front is an 
issue that involves real complications for the three-body 
system. For an extensive discussion of the problem see 
Ref. II^. Presently, we treat the spin in a simplified 
manner; namely by averaging over elementary spin pro- 
jections. This approximation leads to an effective one- 
channel bose-type three-body equation. This simplifies 
the calculation significantly for the sake of a transparent 
treatment of in-medium effects of the relativistic equa- 
tions. A bose-type relativistic three-body equation with 
a zero range interaction has previously proven useful for 
studies of isolated nucleon form factors [ pT|j32| . The tech- 
nically involved full inclusion of the spin degrees of free- 
dom along the lines of Ref. |^ will be postponed to a 
later stage of the calculations. As a consequence we are 
presently not able to distinguish between so-called scalar, 
pseudoscalar, vector, and axial vector quark pairs. As 
many aspects can be improved, we emphasize that we 
present and solve for the first time an in-medium rela- 
tivistic three-body equation for both finite temperatures 
and finite densities. 

The light front has been utilized earlier for infinite nu- 
clear matter calculations at zero temperatures and nor- 
mal nuclear matter densities in ||33|] . That approach is 
based on the Briickner G matrix and two-body correla- 
tions are treated with an in-medium Blankenbecler-Sugar 
equation including a medium-modified interaction. For a 
general discussion of differences between the Green func- 
tion approach (for finite temperatures and densities) used 
here and the Bruckner approach (at zero temperatures) 
see e.g. Ref. Q. 

The equations derived here are based on a systematic 
quantum statistical framework using a cluster expansion 
for the Green functions. We start with the chronological 
Green function MM 



^Qlj/ = {TrA^[r)Al{T')) 

^6{t-t'){A,At)A\{t')) 
T6{T'-T){Al{T')A^{r)) 



(1) 



where t denotes the light cone time = t + z. The up- 
per (lower) sign stands for Fermions (Bosons) , A^ (r) are 
arbitrary Heisenberg operators. For the three-particle 
problem in question they are given by Dirac operators 



taken at equal r, viz. A^ij) = V'i(''')V'2(''')V'3 ('''); indices 
represent all other quantum numbers. Note that on the 
light front the chronological Green function for spin 1/2 
(in this case Aa(r) = V'l (''")) differs from the Feynman 
propagator by a contact term. Presently we neglect this 
term being of higher order |^,^. The generalization 
to imaginary time based on the Matsubara techniques 
is direct and has been used, e.g., to treat few-body cor- 

19)). The 
by § 



relations in nucleonic systems (see e.g. [|9|,p^l 
equation of motion (Dyson equation) is given 
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where 



A/'a/3 = (K(r),4(r)]±) 



(2) 



(3) 



Here, we neglect retardation in the mass operator A^, 
that would introduce more intermediate Fock space com- 
ponents, viz. 

Ml-/ ^ 5{t - T')Ml^p, (4) 
= E([K,i/](r),v4t(r)]±)Ag^ (5) 



So far the formalism presented is rather general and for a 
typical many-body Hamiltonian (with generic two-body 
interactions V2) has been proven useful for various do- 
mains of many-body physics (see ) including the calcu- 
lation of correlations in nuclear systems at finite densities 
and temperatures []rT|-|l9| . 

At chemical potentials /i and temperatures T — 
l/lksP) averaging is due to a grand canonical ensemble 
in equilibrium, (O) = TrjpoO}, where po denotes the cor- 
responding statistical operator expressed in terms of light 
front operators. In averaging we consider the few-body 
cluster embedded in a homogeneous mean field of un- 
correlated particles which is a reasonable approximation. 
Hence, together with (||) the hierarchy of Green functions 
decouples. This leads, upon introducing particle b{k) and 
antiparticle d{k) Fock operators, to the standard Fermi- 
Dirac single particle distribution functions. For particles 
this reads 



/(fc+,fci) = (6t(fc)fe(fc)) 



exp 



[3 { -{k+ + k-,) - 



(6) 



expressed in terms of light front form momenta given by 
k± = {kx,ky) and k^ = ko ± k^. The on-shell light- 
front energy is given by fc"^ — (^1 + "^(M; TY)/k'^ with 
a medium dependent mass m{fi,T). The components of 

the four vector /con are fcon — (fcoiu ^-L' k^)- The thermo- 
dynamic Green function for the one-particle case, by use 
of eqs. (|l|), (H), and (^) is then given by 



2 



1-/ 



2 '^on 



(7) 



Presently we consider only particle degrees of freedom 
on the light front and have therefore dropped the term 
related to 0{—k^). The energy variable k~ plays the 
role of the Matsubara frequency ||], ^{k^ + fc+) = fc° ^ 
7rA/(— + fi. Eventually, averaging over the elementary 
spin projections leads to a bose type Green function. 

For particles that are part of a larger n-body cluster it 
is convenient to introduce fractions x = k~^ / , where 
P^ is the plus component of the cluster's cm. momen- 
tum. For a cluster at rest P+ = M„, where Mn is the 
mass of the cluster. 



FIG. 1. Feynman-Galitzkii equation for the two-body 
t-matrix with zero range interaction. The crosses refer to 
the Pauli-blocking factor Lines represent quasi-particles. 

Evaluation of Eq. (j^) for two or three particles with 
the mass operator of eq. (j^) in a homogeneous medium 
of independent particles leads to resolvent equations of 
the Green function where the single time r — r' is Fourier 
transformed to a Matsubara frequency. These equations 
include mass corrections and Pauli blocking factors in a 
systematic way. They can be rearranged into equations 
for t-matrices. The resulting equation for the two-body t- 
matrix T2 has the same formal structure as the Feynman- 
Galitzkii equation shown in Fig. ^ 



T2^V2 + V2R0N2T2 



(8) 



where V2 represents the two-body interaction, and Rq the 
interaction independent two-body resolvent. The Pauli 
blocking factor 7V2, represented by the crosses in Fig. |^, 
is given by 



N2 = /1/2 - /1/2, f^l-f: 



(9) 



where the indices of the Fcrmi-Dirac function / reflect 
particle quantum numbers. 
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FIG. 2. Loop diagram corresponding to the kernel of the 
integral equation (^|) . The crosses refer to the Pauli-blocking 
factor 



For a numerical analysis and a first calculation of in- 
medium effects we utilize a zero-range model studied ear- 
lier in a different context |^,^. The equation repre- 
sented by Fig. ^ can be summed and leads to a solution 
for the two-particle propagator t{M2), i.e. 



r{M2) = {iX-^-B{M2)y 



(10) 



The expression for B{M2) is represented by the loop di- 
agram of Fig. H and, in the rest system of the two-body 
system F'' = (M2, 0, 0, 0), given by 



B{M2) 
where 



dx(fk±_ 1 - F{x,k\) 



(27r)3 J x{l - x) 



{ki+m^)/x{l~x), 



F{x,kl)^ f{x,kl) + f{l-x,kl) 



(11) 



(12) 
(13) 



and / given in eq. (^ with x — k^/P^. The integral 
can obviously be separated into a blocking independent 
term Bo{M2) and the change AB{M2) depending on the 
Fermi functions F{x, k'j^). The integral involving Bq{M2) 
has a logarithmic divergence that can be absorbed in a 
redefinition of A. The physical information introduced 
in the rcnormalization of the amplitude is the mass of 
the two bound particles, M2b{i^,T). If we assume that 
the two particle amplitude t has a pole for M2(/i,T) = 
M2Bi^J■^T) we may write 



B{M2b) = Ba{M2B) + AB{M2b) 



(14) 



where the dependence on T and fi is suppressed in the no- 
tation. Although, we use eq. (14) for finite temperatures 
and densities, in principle, it is enough to assume such a 
bound state for the isolated case only, and calculate the 
corresponding in-medium bound state M2b{p-,T) from 
the above equations. This would require a definite reg- 
ularization procedure which we presently want to omit. 
The assumption of a bound state implies no restriction 
on the conclusions of our investigation concerning the rel- 
ative importance of two- and three-body correlations at 
a given temperature and density. This will be obvious 
from the results presented and discussed at the end of 
this Letter. 

For M2(P2 = 0,fi,T) < 2fi (for ^ < to) and 2m < 
^2(^2 = 0,^,T) < 2fi (for fj, > m) instabilities of 
the Fermi gas against formation of Bose condensation or 
Cooper pairing can occur. As a consequence the distri- 
bution functions need to be modified to include the gap 
energy. The critical temperature Tc is given by the condi- 
tion M2(P2 = 0, ^, Tc) — 2/i, see e.g. for a discussion 
of similar aspects in nuclear matter. Implementation of 
such effects into the three-body problem would require 
the notion of condensate/pairing at finite P2 which is a 
difficult problem that has hardly been investigated. We 
presently neglect this as well as the appearance of other 
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two-body correlations in the Fermi functions when solv- 
ing the three-body problem. However, in its turn, the 
Dyson approach (cluster mean field expansion) utilized 
here has the potential to re-evaluate the condition of 
critical temperature in the presence of strong three-body 
correlations; that is, however, not the goal of our present 
work. 

The subtraction imposed by condition (|lj) in the de- 
nominator of eq. ( p^ ) makes r(M2) finite. Note that 
AB{M2b) for large momenta k±, so that AB{M2b) 
is finite. The resulting expression for the two-body prop- 
agator is then given by 



r(M2 



(^i [K(A/2s)arctan2K ^{A-'hs) 
-K(M2)arctan2K"^(A/2)] /(27r) 
+AB{M2b) - ABiAh)] 



where 



KiM2) 



Mi 



(15) 



(16) 



The quark mass m(/i, T) as well as the two-body masses 
M2{fJ-,T), Af2B(M,T), and Af2o(M,T) depend on the 
chemical potential and the temperature of the medium. 
Also there is an additional dependence on the cluster's 
cm. momentum related to the momentum dependence of 
the blocking factors. Hence several changes of kinemat- 
ical variables are necessary in the two-body amplitude 
before eq. ( |l5| ) may be employed in a three-body cluster. 
They are discussed in the following. The two-body mass 



Mi = (Pa - 



= (Af3 - q+) Af3 - 



(17) 



where q denotes the momentum of the odd particle. The 
three-body system is taken at rest, = (Afa, 0,0,0). 
Also the blocking dependent part AB needs to be re- 
visited, since the blocking factors depend on the cluster 
momentum and hence M2 is different for a moving sys- 
tem not only because of eq. ( p^ ) , but also due to medium 
effects that depend on the momentum. As a consequence 
the arguments of the blocking factors of eq. (|l^) have to 
be properly replaced by 

F{x, y- k^,qi_) = /(x, k^) + /(I -x-y, {k + q)l). 

(18) 

In the three-body rest frame x — k^ jM-j, and y = /M^. 

The homogenous AGS-type in-medium equations for 
finite temperatures and densities for the nonrelativistic 
bound states have been given in Refs. |l|,|l9|. A dia- 
grammatic representation of these equations for a zero 
range interaction is given in Fig. ||. Based on the Fock 



space representation, a derivation of relativistic three- 
body equations on the light front is formally identical to 
the nonrelativistic case, if we neglect antiparticle degrees 
of freedom. It appears that the modifications required 
to arrive at an in-medium relativistic three-body equa- 
tion for (spin averaged) quasi particles on the light front 
are close to the modifications needed in the nonrelativis- 
tic case. The reason is the formal similarity of the light 
front and nonrelativistic approaches after the approxi- 
mations that have been mentioned above. Hence a rel- 
ativistic AGS-type equation on the light front including 
the effects of finite temperature and density can then be 
written as (compare Ref. [^Ip^ for the isolated case) 



FIG. 3. Diagrammatic representation of the in-medium 
AGS-type equation for a zero range interaction. The crosses 
refer to the Pauli-blocking factor N2- The two-body input is 
given in Fig. ^ and Eq. |lo| 





r(2;,'f±) 



t{M2) 



dx 



(27r)3 ■ JM^Mi x{l~y- x) 
,2 I - F{x,y;k^,q^) 



d^k 



Mi - Mi, 



r{x,k±_), 



(19) 



where we have introduced vertex functions F and t{M2) 
given before. Here 



kT"" = ^{l~x){xMi-m^), 



(20) 



and the mass of the virtual three-particle state (in the 
rest system) is 



Mi 



k\+m? (fi+m? {k + q)\ + m? 



V 



1-x-y 



(21) 



which is the sum of the on-shcU minus-components of the 
three particles. 

The equations suggested here, which result from a sys- 
tematic Dyson equation approach, differ formally from 
those of Ref. Q. For the case T = at finite densi- 
ties the blocking factors could be replaced by /(/i,T) 
n{kp{^)) = 9{k — kp) where in the notation of Ref. 0| 
/i(T = 0) = ^ykp + m^. In Ref. all three quark mo- 
menta are restricted according to eq. (p2|) 



(l-ni)(l-n2)(l-n3) 



(22) 



instead of 1 — ni — n2 (and permutations) from eqs. ( |lq ) 
and (^9|) . Note that in these three-body equations an ad- 
ditional blocking of the spectator particle does not arise 
in a three-body equation driven by a two-body interac- 
tion. 
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FIG. 4. Two-quark vs. three-quark binding energies, 
B„{fM,T)/m{fi,T), at T = 10 MeV for different chemical po- 
tentials ^ as indicated. Dotted line: M2/2 — M3/3. 

Presently, instead of deriving and solving two- and 
three-body equations for one specific model, our inten- 
tion is to explore the correlations between two- and three- 
particle binding energies in a medium of finite temper- 
atures and densities and explore the three-body Mott 
transision. For the time being we chose values of m{^, T) 
provided earlier by a Nambu- Jona-Lasinio model |^,^ . 
That model's approximations are close to those used in 
the present solution of the relativistic three-body equa- 
tion with a simple zero range interaction. The values 
used for this calculation are given in Table 0. Quark 
masses are chosen for values of /x and T to study the 
Mott transition. 



TABLE I. Quark masses m(/i, T) for different tempera- 
tures T and chemical potentials /i used in the calculation |^ . 
All values in MeV. For the isolated case m = 300 MeV. 



/i\r 


10 


50 


100 


100 


300 


299 


279 


200 


300 


293 


236 


300 


289 






307 


279 







TABLE IL Three-body binding energies for weak coupling 
{B2/m = 10"^) and ^ = 100 MeV. 



r[MeV] 10 50 100 

B3[MeV] 23 19 5l) 

Bz/m[%] 7.6 6.5 1.8 



The correlations between the two-body and three-body 
binding energies, B2 and are shown in Fig. ^ for a 
temperature of T = 10 MeV. The binding energies are 
defined by 

B3{fi,T) = m{^i,T) + M2B{f^,T) - M3B{f^,T) (23) 
B2(/i,T) = 2m{p,T) - M2BUhT). (24) 

Binding energies are shown in units of the respective 
quark masses at the temperature and chemical poten- 
tial indicated. The dotted line indicates M2/2 = M3/3 
(compare ) and in a simple chemical picture of an ideal 
gas (law of mass action) the relative importance of the 
clusters. The isolated case (i.e. no medium) is shown 
as a solid line. As the chemical potential increases (and 
the quark mass becomes smaller) three-body correlation 
become weaker for a given B2. In this particular case 
T = 10 MeV the three-body bound state disappears for 
ji ~ 310 MeV, even for a model that allows for a two- 
quark bound state. 

In the weak coupling case B2 ^ Q three-body correla- 
tions appear stronger, however, for a given temperature 
the corresponding bound states vanish above a certain 
value for the chemical potential. For a given chemical 
potential /i = 100 MeV and in the weak coupling limit 
the respective values for B^, are given in Table |l|. 

The dependence of B3 (T) for a given chemical poten- 
tial is opposite from the nuclear case, compare Ref. [Q. 
This is in agreement with the expectation from other ap- 
proaches that give a negative slope for the phase transi- 
tion, see e.g. For strong couplings, i.e. B2/m > 0.65 
and T = 10 MeV, two-body correlations dominate for all 
chemical potentials. 




B/m 

FIG. 5. Two-quark vs. three-quark binding energies for 
different chemical potentials fi and temperatures T = 100 
MeV (dashed hues) and T = 50 MeV (dashed-dotted hues). 
Dotted line: M2/2 = M3/3. 
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In Fig. H the same correlations are shown for different 
temperatures and chemical potentials. For higher tem- 
peratures and chemical potentials the three-body bound 
state disappears and also the three-body correlations be- 
come weaker as the lines are below the long-dashed line. 
In all examples the isolated case is limiting. An increase 
in fi means more particles in the medium that prefer- 
entially occupy the momentum components necessary to 
form bound states. A similar effect arises for increasing 
temperatures. For a given chemical potential = 100 
MeV Fig. ^ reflects stronger correlations at lower tem- 
peratures. 
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FIG. 6. Binding energy of the three-body system for 
M2/m = 1 for different temperatures as indicated. B3 — 
corresponds to the Mott density. 

In Fig. H we show the dependence of the three-body 
binding energy on the chemical potential for different 
temperatures T for an assumed two-body bound state 
of M2/m = 1. The continuum (for 2q+lq break-up) is 
reached at B3 = for certain chemical potentials known 
as Mott transition. For increasing temperature the chem- 
ical potential where the Mott transition occurs becomes 
smaller. 

The Mott lines, i.e. the values of T and /i where 
the transitions occur are given in Fig. |^ which is our 
main result. This is in qualitative agreement with the 
confinement-deconfinement phase transition |^ . In these 
cases we have assumed a value for the two-quark bound 
state of M2/TO = 1, which is close to a value predicted 
by the NJL model |^ and for a lesser bound system of 
M2/TO = 7/4 to study the sensitivity on the two-body 
input. It is obvious that the qualitative behavior of the 
Mott transition retains for the two different models. 

The onset of superfluidity is expected for M2{P2 — 
0,fj.,Tc) — 2/i. We mention here that for e.g. for Tc — 



10 MeV this condition is fulfilled at /i = 150 MeV for 
M2 = m and at = 260 MeV for the case M2 = 7m/ 4. 
However, no definite conclusion can be made on the basis 
of the present treatment of the two-body amplitude. 



200 




400 



FIG. 7. Mott line for the three-body system at rest in the 
medium with M^/m — 1 (solid) andM2/m = 7/4 (dashed). 
For values of T and fi below the Mott lines three-body bound 
states can be formed. 

In conclusion, we have derived for the first time a con- 
sistent relativistic three-body equation for particles em- 
bedded in a medium of both finite temperature and finite 
density. This equation systematically includes the effects 
of Pauli blocking and mass shift of the quasi-particles in- 
volved. The equations are solved for parameters close to 
the phase transition of QCD. We find that three-body 
clusters become less stable for a denser or hotter system 
than the two-body cluster. This justifies, a posteriori 
our restriction to the two-body bound states in this in- 
vestigation. However correlations may still exist as the 
"pole" moves into the continuum (e.g. corresponding to 
anti-bound states in this approximation). 

For further investigations and to achieve more quanti- 
tative results a specific model can be chosen that leads to 
a specific two-body t-matrix and a corresponding three- 
body bound state. Furthermore, the approximations we 
used in the treatment of spin of the particles have to be 
relaxed and a full treatment needs to be implemented 
along the lines discussed at length in Ref. |^^. The ne- 
glect of confinement in the vicinity of the Mott line may 
be justified by color screening discussed in Ref. [Q. How- 
ever a more realistic treatment of confinement is highly 
desirable . The light front approach leading to 3-dim 
equations is suited to include relativistic confining poten- 
tials into the formalism. 

In view of the present simplifications, related to the av- 
eraging of spin degrees of freedom and the use of a zero 
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range force, we would like to emphasize that the main 
result of this Letter is to show that an in-medium three- 
body equation that reflects the basic requirements of rel- 
ativity can be derived and solved in the physical region 
of interest close to the Mott transition. The light-front 
framework appears to be very useful in this context as it 
allows for a Hamiltonian formalism close to an intuitive 
interpretation based on nonrelativistic approaches. 

Addressing the question of color-superconductivity in 
a next step the full three-body t-matrix has to be im- 
plemented in a calculation of the critical temperature of 
condensation/pairing. This can also be achieved in a sys- 
tematic fashion within the Dyson equation approach used 
here. 
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